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It was recently proposed that the stress state of a material can also be altered via electron or hole doping, 
a concept termed electronic stress (ES), which is different from the traditional mechanical stress (MS) due to 
lattice contraction or expansion. Here we demonstrate the equivalence of ES and MS in structural stabilization, 
using In wires on Si(lll) as a prototypical example. Our systematic density-functional theory calculations 
reveal that, first, for the same degrees of carrier doping into the In wires, the ES of the high-temperature metallic 
4x1 structure is only slightly compressive, while that of the low-temperature insulating 8x2 structure is much 
larger and highly anisotropic. As a consequence, the intrinsic energy difference between the two phases is 
significantly reduced towards electronically phase-separated ground states. Our calculations further demonstrate 
quantitatively that such intriguing phase tunabilities can be achieved equivalently via lattice-contraction induced 
MS in the absence of charge doping. We also validate the equivalence through our detailed scanning tunneling 
microscopy experiments. The present findings have important implications in understanding the underlying 
driving forces involved in various phase transitions of simple and complex systems alike. 

PACS numbers: 73.20.At, 68.35.Md, 71.30.+h 


I. INTRODUCTION 


Mechanical stress (MS) produced by lattice deformation 
is well established to tune the electronic, magnetic, optical, 
and phononic properties of material^, and such an elastic 
stress (strain) engineering has been widely adopted to sub¬ 
stantially improve the carrier mobilities in semiconductor de¬ 
vices^. These MS-driven tuning effects are more profound in 
strongly-correlated or low-dimensional materials, mainly due 
to the enhanced entanglement between lattice, charge, spin, 
and orbit degrees of freedom. For example, MS has been 
demonstrated to tune the Mott transition temperature in VO 2 
nanowires^, and to generate giant pseudomagnetic field and 
band gap in graphene^. The creation of MS usually relies on 
high-pressure instruments or lattice mismatch engineering at 
the interfaces. Contrasting with the MS induced by lattice de¬ 
formation, the so-called quantum electronic stress (ES), a pure 
electronic effect on the stress originating from the variation of 
carrier density, has been recently introduced and formulated 
within density functional theory (DFT)&. Indeed, the ES in¬ 
duced by quantum electronic confinement in metal thin films 
has been demonstrated theoretically^ - - and experimentally^. 
Since the MS and ES have substantially different origins in¬ 
volving explicitly the variations of lattice and charge degrees 
of freedom, respectively, it is interesting and challenging to 
examine whether and how they can equivalently tune the phys¬ 
ical properties, especially in the same system. 

For the surface structures formed by epitaxial metal atom 


adsorption on semiconductor surfaces, there are frequently 
competing electronic phases^ 13 because of their reduced 
phase space, and the stability of these phases can be effec¬ 
tively tuned by deforming the lattice^ or by varying charge 
carriers^. In this sense, the low-dimensional electronic sys¬ 
tems formed on surfaces provide a unique playground to 
demonstrate the tuning effect of phase stability in terms of 
surface MS and ES (hereafter MS and ES refer to the surface 
ones). Here, we focus on a prototypical example of quasi-one 
dimensional (ID) systems, self-assembled Indium (In) atom 
wires on the Si(lll) surface (see Fig. 1). This In/Si(l 11) 
surface system undergoes a structural phase transition from 
a high-temperature metallic 4x1 phase [see Fig. 1(a)] to a 
low-temperature insulating 8x2 phase [Fig. 1(b)] at a tran¬ 
sition temperature ( T c ) of —125 K—. The structural model 
of the 8x2 phase is well established to have the formation 
of In hexagons via a periodic lattice distortion^ - - 7 , reflecting 
the presence of the MS. To produce and quantify the ES that 
may influence the relative stability of the 4x1 and 8x2 struc¬ 
tures, we introduce electron doping into In wires. It is very in¬ 
teresting to explore how such an electron-doping induced ES 
changes depending on the metallic and insulating phases, and 
also to examine the equivalence of the ES and MS in tuning 
the phase stabilization of the 4x 1 and 8x2 structures. 

In this paper, we present a comprehensive study of surface 
MS and ES in the In/Si(l 11) system by using a van der Waals 
(vdW) energy-corrected hybrid DFT calculation. We find that 
the formation of In hexagons brings a significant reduction of 
the tensile MS perpendicular to In wires, leading to the sta- 
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FIG. 1: (Color on line) Top view of the optimized (a) 4x1 and (b) 
8x2 structures of In/Si(l 11). Each In atom wire is composed of two 
zigzag chains of In atoms. The dark and gray circles represent In and 
Si atoms, respectively. For distinction, Si atoms in the subsurface are 
drawn with small circles. Each unit cell is indicated by the dotted 
line. The x (y) axis is parallel (perpendicular) to In chains. Numbers 
denote the In-In distance (in A) between In chains. The arrows in (b) 
indicate schematically the relaxation directions of In atoms forming 
the In hexagon (see Table I), when electrons are doped to In wires. 


bilization of the 8x2 structure. Interestingly, the ES induced 
by electron doping into In wires exhibits drastically different 
features between the 4x1 and 8x2 structures: i.e., the ES of 
the 4x1 metallic structure is slightly compressive, while that 
of the 8x2 insulating structure is anisotropic with a highly 
compressive (tensile) component along the direction parallel 
(perpendicular) to In wires. As a result, the surface energy dif¬ 
ference between the 4 x 1 and 8x2 structures decreases as the 
amount of electron doping increases. This ES-driven tuning 
effect on the relative stability of the 4x1 and 8x2 structures 
is found to be equivalent to the MS-driven one obtained by 
applying a compressive lattice strain. Our theoretical predic¬ 
tions are confirmed by scanning tunneling microscopy (STM) 
measurement at 5 K. 


II. CALCULATIONAL METHOD 


The present hybrid DFT+vdW calculation^^ was per¬ 
formed using the FHI-aims^i code for an accurate, all-electron 
description based on numeric atom-centered orbitals, with 
“tight” computational settings. For the exchange-correlation 
energy, we employed the screened hybrid functional of Heyd- 
Scuseria-Ernzerhof (HSE) 21 i 22 . The Si(lll) substrate (with 
the Si lattice constant ao = 5.418 A) was modeled by a 6- 
layer slab (not including the Si surface chain bonded to the In 
chains) with ^30 A of vacuum in between the slabs, where 
each Si atom in the bottom layer was passivated by one H 
atom. The k-space integrations were done equivalently with 
64 and 16 k points in the surface Brillouin zone of the 4x 1 
and 8x2 unit cells, respectively. All atoms except the bottom 
layer were allowed to relax along the calculated forces until 
all the residual force components were less than 0.001 eV/A. 
The employed hybrid DFT+vdW scheme was successfully ap¬ 
plied not only to determine the energy stability of the 4 x 1 and 
8x2 structures of In/Si(l 11) 17 but also to calculate the stress 
tensor—. 


III. RESULTS 


We begin to optimize the 4x1 and 8x2 structures without 
electron doping by using the hybrid DFT+vdW scheme. The 
optimized 4x1 and 8x2 structures are displayed in Fig. 1(a) 
and 1(b), respectively. It is seen that the 8x2 structure has the 
shorter In-In distances (<7i n -in = 2.957 and 2.962 A) between 
two In chains compared to that (3.085 A) in the 4x1 struc¬ 
ture, forming In hexagons. Such an 8x2 hexagon structure is 
found to be more stable than the 4 x 1 structure by 33 meV per 
4x1 unit cell. The calculated surface band structures of the 
4x1 and 8x2 structures show that the 4x1 structure exhibits 
the presence of three metallic bands crossing the Fermi level 
whereas the 8x2 structure has a band gap of 0.31 eV, in good 
agreement with previous experimental data^i - —. To examine 
how the MS changes after the formation of In hexagons, we 
calculate the MS difference Ag^ between the 8x2 (a) and 
4x1 (P) structures, defined as 


A M M M 
AG ij ~ C ij-u G ijfi 


1 3(A a y a ) 1 3(Apyp) 

Aa d&i j Ap d&ij 
1 dE s iab : a 1 
Aa dZi j Ap dZi j 


( 1 ) 


Here, Zij (/, j = x,y) denotes the element of strain tensor, A the 
surface area, and y (E s i a b) the surface (slab) energy. For the 
last equality, see the Appendix. Therefore, Acan be eval¬ 
uated by using E s i a b obtained from the slab calculation. The 
calculated results of Aare plotted in Fig. 2(a), together 
with those (discussed below) obtained with electron doping. 
We find that (i) the G^[ and G^ components in 8x2 are al¬ 
most the same as those in 4 x 1 and (ii) the G^ component in 
8x2 is reduced as much as 29.79 meV/A 2 compared to that 
in 4x 1. Thus, the In-hexagon formation results in a decrease 
in the tensile surface stress along the y direction, giving rise 
to the stabilization of the 8x2 structure. Here, we note that 
the absolute value of G^ in the 4x 1 reference is G^ = 54.79, 
G^ = 121.68, and G^ = 0 meV/A 2 , indicating a tensile surface 
stress along the x and y directions. 

Next, we study the ES induced by electron doping with ex¬ 
cess electronic charge n e per 4x1 unit cell 27 . The concept of 
ES was recently formulated within DFT-, and it can be practi¬ 
cally calculated by using the difference of the MS obtained at 
the total electronic charge = no + n e and that at the ground- 
state electronic charge no : 


Gfj(n e )=^(n d )-^(n 0 ) 


1 3 E s i a b{nd ) 1 dE slab (no) 

A dZij A dz ij 


We here consider the two different states for the treatment of 
G fj( n d) : one is the “initial” state without the relaxation of 
atoms (i.e., fixing the structure having no) and the other is the 
“final” state which allows the atomic relaxation along the gen¬ 
erated forces due to electron doping. The calculated initial- 
state and final-state ES results for the 4 x 1 and 8x2 structures 
are plotted as a function of n e in Fig. 2(b) and 2(c), respec¬ 
tively. For the 4x1 initial state, G^ x (g^) is slightly nega¬ 
tive as -1.15 (-0.01), -1.15 (-1.81), and -1.16 (-3.96) 
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FIG. 2: (Color on line) (a) Calculated MS difference A ofj between the 8x2 and 4x1 structures as a function of n e . The ES afj for the 4x 1 
and 8x2 structures are given in (b) and (c), respectively. 


meV/A 2 for n e = 0.05, 0.1, and 0.15e, respectively. The in¬ 
clusion of lattice relaxation within the 4x1 final state shows 
a negligible change in ofj [see Fig. 2(b)]. Thus, we can say 
that electron doping in the 4 x 1 structure produces a weakly 
compressive ES. Contrasting with the 4x1 case, the 8x2 ini¬ 
tial state exhibits larger electron-doping effects with (o^) 
= -3.50 (+4.29), -8.34 (+6.38), and -9.22 (+6.39) meV/A 2 
for n e = 0.05, 0.1, and 0.15^, respectively. As shown in Fig. 
2(c), the 8 x2 final state further increases the magnitude of o xx 
(G^) as -4.07 (+7.91), -9.55 (+14.30), and -10.34 (+17.25) 
meV/A 2 for n e = 0.05, 0.1, and 0.15e, respectively, and their 
magnitudes monotonically increase with increasing n e . Not¬ 
ing that the In-hexagon formation in the 8x2 structure in¬ 
volves a reduction of the mechanical tensile surface stress [see 
Fig. 2(a)], the significant final-state effect of a yy may accom¬ 
pany a large atomic relaxation. Indeed, Table I shows that 
the electron doping of n e = O.le in the 8x2 structure gives a 
conspicuous relaxation of In atoms forming the In hexagon, as 
indicated by the arrows in Fig. 1(b). On the other hand, for the 
4x1 structure, there is a negligible atomic relaxation caused 


TABLE I: Calculated displacements (in A) of In and Si atoms in the 
4x1 and 8x2 final-state structures with n e = O.le relative to the po¬ 
sitions obtained from the corresponding structures without electron 
doping. The labeling of In and Si atoms is shown in Fig. 1. 



Ax 

4x1 

Ay 

Az 

Ax 

8x2 

Ay 

Az 

Ini 

0.000 

0.001 

-0.007 

0.043 

0.000 

-0.010 

In 2 

0.000 

-0.002 

0.007 

0.009 

0.071 

0.022 

I n 3 

0.000 

-0.006 

0.009 

-0.011 

0.046 

-0.025 

In4 

0.000 

-0.008 

-0.014 

0.132 

-0.011 

0.003 

I n 5 




-0.124 

0.017 

0.008 

In 6 




0.010 

-0.044 

-0.032 

In 7 




-0.008 

-0.062 

0.024 

In 8 




-0.044 

0.006 

-0.006 

Sii 

0.000 

-0.004 

-0.001 

0.010 

-0.005 

-0.002 

Si 2 

0.000 

-0.003 

-0.004 

-0.006 

0.015 

0.010 

Si 3 




-0.008 

-0.015 

0.006 

Si4 




0.003 

0.014 

0.001 


by electron doping (see Table I). It is remarkable that the ES 
of the 8 x 2 structure is anisotropic with a highly compressive 
(tensile) component along the direction parallel (perpendic¬ 
ular) to In wires, and thus their magnitudes are significantly 
larger than those of the 4x 1 structure. 

To account for the different features of ES between the 4 x 1 
and 8x2 structures, we display in Fig. 3 the charge charac¬ 
ters of their final states with n e = O.le, together with the cor¬ 
responding band structures. For the 4x1 structure, electron 
doping shifts the Fermi level continuously upward by occu¬ 
pying the metallic states, and the occupied excess electrons 
are found to be well distributed over the whole In wires [see 
Fig. 3(a)], indicating a widely delocalized metallic character. 
On the other hand, for the 8x2 structure, the excess electrons 
occupying the conduction bands above the band gap show a 
strongly delocalized character along each chain with charge 
depletion between In chains [see Fig. 3(b)], leading to en¬ 
hance the metallic bonding along each chain while to weaken 



FIG. 3: (Color on line) Calculated surface band structures of the 
electron-doped (a) 4x1 and (b) 8x2 final states with n e = O.le per 
4x1 unit cell. The energy zero represents the Fermi level £>. The 
surface Brilloiun zone is displayed in (b). For each structure, the 
charge character of excess electrons, obtained by the charge density 
difference p Hd — p„ 0 , is also given. Here, the charge contour plot with 
a contour spacing of 0.2 x 10 _3 e/A 3 is drawn in a lateral plane near 
In atom wires. 
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FIG. 4: (Color on line) Calculated total-energy difference 
AFg X 2-4x l per 4x 1 unit cell between the 8x2 and 4x1 structures as 
a function of electron doping n e . The arrow indicates the heavy elec¬ 
tron doping of ^O.le in the STM and angle-resolved photoemission 
spectroscopy. 


the strength of covalently bound In hexagons formed by In 2 - 
L 15 and In 4 -In 7 covalent bondsii Based on such contrasting 
charge characters of excess electrons between the 4x1 and 
8x2 structures, it is likely that the 4 x 1 structure has a weakly 
compressive ES, while the 8x2 structure has an anisotropic 
feature of ES with a highly compressive (tensile) stress along 
the x (y) direction. Here, the electron-doped 8x2 structure, 
which weakens the covalent bonding of In atoms between In 
chains, in turn gives an increase of tensile stress along the y 
direction. 

By summation of G^(no) and cfj(n e ), we can obtain the 
MS at rid - no + n e . Accordingly, the difference of 

MS between the electron-doped 8x2 and 4x1 structures is 
given by (rid) = Ag-^ (no) + A afj(n e ). The calculated re¬ 
sults for A o^j(nd) are plotted as a function of n e in Fig. 2(a). It 
is seen that the slope of increase in Aa^ is greater than that of 
decrease in Ag^J, thereby giving rise to a decrease in the mag¬ 
nitude of Ag^+Ag^ with increasing n e . Consequently, one 
expects a reduction of the total-energy difference A^sx 2 - 4 x 1 
(per 4x1 unit cell) between the 8 x 2 and 4x1 structures with 
increasing n e . Indeed, as shown in Fig. 4, the present calcula¬ 
tion of AEg x 2 - 4 x 1 shows that the relative stability of 8 x 2 with 
respect to 4x 1 decreases with increasing n e . Interestingly, the 
4x1 structure becomes more stabilized than the 8 x 2 struc¬ 
ture above n e ~ O.lle, implying that the ground state of the 
In/Si(l 11) system can be switched to the 4x1 structure by 
electron doping. 

It is noteworthy that the decrease in the magnitude of 
AEg x 2—4x1 with electron doping is consistent with several 
previous experimental observations that (i) ft-type Si(lll) 
substrate yields the coexistence of the 4 x 1 and 8x2 phases 
even at 47 K, whose areal ratio can be tuned by optical excita¬ 
tion that decreases the amount of electron doping in In wiresS 
and (ii) electron doping via Na adsorption on the In/Si(l 11) 
surface suppresses the (4x l)«-»( 8 x 2 ) phase transition, result¬ 
ing in a lowering of T c —. Here, electron doping with the 
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FIG. 5: (Color on line) (a) Calculated MS components o^J, o^, and 
for the 4x 1 and 8x2 structures and (b) MS difference AafJ be¬ 
tween the 8x2 and 4x1 structures as a function of \Aa\/aQ. 


Na coverage of ~0.02 ML preserved the 4x1 phase even at 
50 K. Remarkably, this Na coverage was estimated to give 
n e & O.l^—, at which our calculated value of A£g X 2 - 4 xi ap¬ 
proaches zero (see Fig. 4)22. 

For comparison with the features of ES, we further study 
the MS induced by lattice deformation. According to our re¬ 
cent STM and DFT studyA2, the vacancy creation in In wires 
produces a compressive lattice strain to change the relative 
stability of the 4x1 and 8x2 structures, as discussed below. 
We here examine how the MS evolves with contracting the 
lattice constant a of the Si(lll) substrate by 1%. As shown 
in Fig. 5(a), we find that the tensile MS components g^J 
and G^J in the 4x1 (8x2) structures decrease with contract¬ 
ing a and are finally converted to be compressive at a con¬ 
traction of ~0.3(0.3) and ~0.6(0.5)%, respectively. Interest¬ 
ingly, the MS difference Ag^ (Ag^) between the 8x2 and 
4x1 structures decreases (increases) with contracting a [see 
Fig. 5(b)], similar to the pattern of Ag^[ (Ag^) as a func¬ 
tion of n e in the above-mentioned case of electron doping [see 
Fig. 2(a)]. The resulting magnitude of Ag^+Ag^ overall de¬ 
creases with contracting a , which in turn decreases the magni¬ 
tude of AE'gx2—4xi as 28.5, 24.1, and 18.1 meV per 4x1 unit 
cell at \Aa\/ao = 0.3, 0 . 6 , and 1 %, respectively. Thus, we can 
say that both the ES induced by electron doping and the MS 
induced by lattice contraction equally contribute to tune the 
structural phase stabilization in the In/Si(l 11) system. 

In order to verify our theoretical prediction of the MS- and 
ES-driven tuning effects on the stabilities of competing 4x1 
and 8x2 phases, we have performed STM experiments at 5 
Ks 2I. We intentionally created vacancy defects in In wires to 
induce strain fields. It is observed that only the 8x2 phase 
exists at low defect density [Fig. 6 (a)], while both the 4x1 
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FIG. 6: (Color on line) STM images of the In atom wires on Si(l 11) 
substrates, (a) Defect-poor In wires on n-type Si (defect density of 
0.0038 nm -2 ), V s = — 1 V, I t = 5 pA. (b) Defect-rich In wires on 
p-type Si (defect density of 0.0085 nm -2 ), V s = — 1 V, I t = 2 pA. 
(c),(d) The same area of defect-rich In wires on n-type Si (defect 
density of 0.0086 nm -2 ), with V s = — 1 Y and 1 V, respectively and 
the same I t of 2 pA. V s and I t denote the sample bias and tunneling 
current, respectively. All the STM images are of the same size and 
were acquired at 5 K. 


and 8 x2 phases coexist at high defect density [Fig. 6(c)]. The 
latter electronically phase-separated ground state can be at¬ 
tributed to large compressive strains^ due to high defect den¬ 
sity, consistent with the present theoretical prediction that the 
magnitudes of Ag^+Ao^ and A£sx 2 - 4 xi decrease with con¬ 
tracting a. On the other hand, the ES-driven tuning effect is 
demonstrated by adopting either n- or p -type substrate: i.e., 
for a certain defect density, only the 8x2 phase is present 
on /7-type substrate (hole doping) [Fig. 6(b)], whereas both 
the 4x1 and 8x2 phases coexist on n -type substrate (elec¬ 
tron doping) [Fig. 6(c)]. Alternatively, we utilize the sur¬ 
face charging effect at 5 K, where carrier relaxation between 
the surface layer and the bulk is substantially suppressed. As 
shown in Fig. 6(c) and 6(d), negative bias voltage (hole dop¬ 
ing) tends to stabilize the 8 x 2 phase, while positive bias volt¬ 
age (electron doping) favors the 4x1 phase. Therefore, it is 
demonstrated that increasing electron doping can favor the 
stabilization of the 4x 1 structure. 


IV. SUMMARY 


We have demonstrated the equivalent roles of ES and MS in 
tuning the relative stability of the 4x1 and 8x2 structures in 
the In/Si(l 11) surface. By means of hybrid DFT+vdW calcu¬ 
lation, we found that electron doping into In wires for the 4 x 1 
and 8x2 structures induces the ES with drastically different 
features, leading to a decrease in the surface energy difference 
between the two structures. We also found that applying a 
compressive lattice strain yields similar results for the surface- 
stress and surface-energy differences between the 4x1 and 


8x2 structures. The equivalent control of phase stability by 
ES and MS has also been validated by low-temperature STM 
experiments. The present findings have important implica¬ 
tions in understanding the underlying driving forces involved 
in various phase transitions of simple and complex systems 
alike, as well as in tailoring the physical properties of such 
systems. 
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APPENDIX: EVALUATION OF THE SURFACE MS 
DIFFERENCE BETWEEN THE 8x2 AND 4x1 STRUC¬ 
TURES 

The surface energy (y) per unit area of the In/Si( 111) sur¬ 
face system can be defined using a centrosymmetric slab ge¬ 
ometry of which both sides consist of two equivalent sur¬ 
faces^: 


Y = ^ IKZ ~ HSiNst - HmNm] , (3) 

where A is the surface area of the unit cell, psi (Min) is the Si 
(In) chemical potential, i.e. the energy per atom in bulk, Nsi 
( Ni n ) is the number of Si (In) atoms in the unit cell, and E s J™ b is 
the total energy of the centrosymmetric slab. The factor of 1/2 
is introduced to take into account the presence of two surfaces 
in the symmetric slab. Thus, the surface energies of the 8x2 
(a) and 4x 1 (P) structures are given as: 


Ya — ^dab,a tysiNsi 4 pi n Ni n \ (4) 

Yp = ^ K,p - WiNsi ~ mNin] • (5) 

The MS difference A off between the 8x2 and 4x1 struc- 
tures can be defined^ as 


A m_ m _ m _ 
AG ij ~ C ij ,p - 


1 3(A a y a ) 

Aa j 


i 

Ap j 


1 as,. 1 agp 


( 6 ) 

( 7 ) 
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where 8 ^ is the surface strain tensor (ij = x,y). The last 
equality holds since the stress of the bulk at the equilibrium 
lattice constant is zero. Assuming that the surface geometry 
obtained using the symmetric slab is the same as that obtained 
using the H-terminated slab, \{E s J^ ha - p) is equal to 

E s iab,a ~ j£E s i ab fi, where E sluh is the total energy obtained 
using the H-terminated slab. Therefore, Eq. (7) can be ex¬ 
pressed in terms of the first order change in the total energy of 


the H-terminated slab: 


1 dE S l a ]j a l dE s l a Jy ^ 

ij - A a d£ij Ap dEij 


( 8 ) 
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